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1. Introduction 

After analysing particle multiplicities in heavy-ion collisions for two decades a remarkably 
simple picture has emerged for the chemical freeze-out parameters [Q, ^, |3]]. Despite much initial 
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Figure 1: Values of the freeze-out parameters obtained at beam energies ranging from 1 GeV to 200 GeV 



skepticism, the thermal model has emerged as a reliable guide for particle multiplicities in heavy 
ion collisions at all collision energies. Some of the results, including analyses from fl EL o, [7fl, 
are summarised in Fig. HI Most of the points in Fig. [I] (except obviously the ones at RHIC) refer 
to integrated (4jt) yields. A clear discrepancy exists in the lower AGS beam energy region be- 
tween the chemical parameters extracted from (published) mid-rapidity yields and those extracted 
using estimates of the 4n yields. The latter tend to give higher values for the chemical freeze-out 
temperature. This will have to be resolved by future experiments at e.g. NICA and FAIR. When 
the temperature and baryon chemical potential are translated to net baryon and energy densities, a 
different, but equivalent, picture emerges shown in Fig. |[ This clearly shows the importance in 
going to the beam energy region of around 8-12 GeV as this corresponds to the highest freeze-out 
baryonic density and to a rapid change in thermodynamic parameters [j|, ^]. 

The dependence of [ig on the invariant beam energy, Js?jn, can be parameterized as ^\ 



Hb(^snn) 



1.308 GeV 



1 + 0.273 GeV 



-l 



'SNN 



Similar dependences have been obtained by other groups [JXJ, Ll|]. and are consistent with the above. 
This predicts that at the LHC pL B « 1 MeV. 

To analyze the changes around 10 GeV use can be made of the entropy density, s, divided 
by T which has been shown to reproduce the freeze-out curve [^] very well. This allows for a 
separation into baryonic and mesonic components, shown in Fig. L3J, it can be seen that mesons 
dominate the chemical freeze-out from about -Jstm rj 10 GeV onwards. 
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Figure 2: The hadronic freeze-out line in the pg — £* phase plane as obtained from the values of jj,b and T 
that have been extracted from the experimental data in [gj . The calculation employs values of jitg and jis 
that ensure (S) = and (Q) = 0.4(B) for each value of /Xg. Also indicated are the beam energies (in GeV/N) 
for which the particular freeze-out conditions are expected at either RHIC or FAIR or NICA. 
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2. Antimatter Production 

One of the striking features of particle production at high energies is the near equal abundance 



of matter and antimatter in the central rapidity region [11, 12]. As is well known a similar symme- 
try existed in the initial stage of the universe and it still remains a mystery as to how this got lost in 
the evolution of the universe reaching a stage with no visible amounts of antimatter being present. 
Closely related to this matter/antimatter symmetry is the production of light antinuclei, hypernuclei 
and antihypernuclei at high energies. Since the first observation of hypernuclei in 1952 [ |I3| ] there 
has been a steady interest in searching for new hypernuclei and exploring the hyperon-nucleon 



interaction which is relevant (see e.g. [ ]14] , |15[ ]) for nuclear physics. Hypernuclei decay with life- 
time which depends on the strength of the hyperon-nucleon interaction. While several hypernuclei 
have been discovered since the first observations in 1952, no antihypernucleus has ever been ob- 
served until the recent discovery of the antihypertriton in Au+Au collisions at Js^n - 200 GeV by 



the STAR collaboration at RHIC [16]. The yield of (anti)hypernuclei measured by STAR is very 
large, in particular they seem to be produced with a similar yield as other (anti)nuclei, in partic- 
ular (anti)helium-3. This abundance is much higher than measured for hypernuclei and nuclei at 
lower energies [[n|]. It is of interest to understand the nature of this enhancement, and for this the 
mechanism of production of (anti)hypernuclei should be investigated. 

The analysis of particle production assessing the degree of thermalization of the particle source 



has been undertaken for many decades [|18|, |19, [20J, 21|, [22]. It has been found that the thermaliza 



tion assumption applies successfully to hadrons produced in a large number of particle and nuclear 
reactions at different energies [23, |4|, 25]. This fact allows us to estimate thermal parameters 



characterizing the particle source for each colliding system, relevant for the understanding of the 
thermal properties of dense and hot matter, and in particular for studies of QCD phase transitions. 
In this paper, using the parametrizations of thermal parameters estimated by the model THER- 



MUS [26, 27] that were shown to best fit the existing data from particle and nuclear collisions 
at several energies, we make thermal model estimates of (anti)hypernuclei that can be directly 
compared to the recently measured unexpected high (anti)hypernuclei yields at RHIC as well as 
predictions of (anti)matter and (anti)hypernuclei production at the Large Hadron Collider (LHC). A 



similar analysis, not including p-p results, has been presented recently in Q28Q where it was shown 
that ratios of hypernuclei to nuclei show an energy dependence similar to the K + /n + one with a 
clear maximum at lower energies. In this paper we study quantitatively how the matter/antimatter 
symmetry is reached as the beam energy is increased. We also estimate ratios of hypernuclei and 
antihypernuclei yields in Au+Au collisions at RHIC using the above mentioned parametrizations 
of thermal parameters that best fit hadron production at RHIC. The present analysis uses a ther- 
mal model and aims to elucidate the production mechanism of hypernuclei and antihypernuclei in 
heavy ion collisions at RHIC and LHC energies, thus providing insight in the surprising increase 
of (anti)hypernuclei production at high energies. 

3. The THERMUS model 

The thermal model assumes that at freeze-out all hadrons in the hadron gas resulting from 
a high energy collision follow equilibrium distributions. The conditions at chemical freeze-out 
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(when inelastic collisions cease) are given by the hadron abundances, while the particle spectra of- 
fer insight into the conditions at thermal freeze-out (when elastic collisions cease). Once evaluated 
the hadron gas partition function gives all primordial thermodynamic quantities of the system by 
simple differentiation. The exact form of the partition function, however, depends on the statistical 
ensemble under consideration. 
Within the grand-canonical ensemble the quantum numbers of the system are conserved on average 



through the action of chemical potentials []24[]. In other words, the baryon content B, strangeness 
content S and charge content Q are fixed on average by }Xb, jJ-s an d \Iq respectively. For each of 
these chemical potentials one can write a corresponding fugacity using the standard prescription 
X = e^' T , where T is the temperature of the system. 

As an example, the density of hadron species i with quantum numbers 5,, S; and Qj, spin- 
isospin degeneracy factor, g/, and mass, m,, emitted directly from the fireball at temperature T is 
given by a second order modified Bessel function of the second kind, 

ni(T,LiB,H S ,HQ, Ys) = ^m^TX*'X s s 'X%Y^K 2 (^). (3.1) 

in the Boltzmann approximation. 

The quantum-statistical result requires either an infinite summation over such K 2 functions or else 



a numerical integration [|26| I27Q . 



The chemical potentials Us and \Lq are typically constrained in applications of the model by 
the initial strangeness and baryon-to-charge ratio in the system under consideration. 

4. Production of antibaryons 

In heavy-ion collisions the increase in the antimatter to matter ratio with the center-of-mass 
energy of the system has been observed earlier by the NA49 [29, ^(J and the STAR [|3l]] collab- 



orations. The trend of p/p ratio increase with the energy towards unity is shown in Fig. |], where 
the open squares refer to heavy ion collisions and the solid circles refer to p-p collisions. It include 
results from the NA49 [||], STAR [^TJ] and the new results from the ALICE Collaboration [JTJ]. 
The resulting baryon chemical potential /i# is shown in Fig. |5] where the dashed line refers to the 



heavy ion description using the THERMUS model [Eq, 27]. The two input parameters, the chem 



ical freeze-out temperature T and the baryon chemical potential /ig as a function of y/s are taken 



fromRef. [32] 



T{ji B )=a-biil-ci4 (4-1) 

with a = 0.166 ±0.002 GeV, b = 0.139 ±0.016 GeV ' and c = 0.053 ±0.021 GeV 3 . This 



parametrization is similar and consistent with the one proposed in Ref. [34]. The solid line in 



Fig. ^ is obtained from THERMUS model [ 26| , 27 ] using T from equation 1 and Hb from equation 
2. The solid circles represent /!#, obtained after fitting experimental data with the THERMUS 



model H2q, 27]. The solid line is a new parametrization adjusted for pp collisions. In view of the 



fact that peripheral and central collisions show no noticeable change in the temperature we have 
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Figure 4: The p/p ratio as function of ^/s. The solid circles are results from p-p collisions and the open 
squares are results from HI collisions as a function of the invariant beam energy p^, p\\ [L4 [30[ pip . 



used the same T dependence for p-p as in heavy ion collisions but the dependence on ji B on beam 
energy is now given by 

HB = d/(l + ey/s) (4.2) 

with d = 0.4 GeV and e = 0.1599 GeV" 1 . 

It is important to note that /»# is always lower in pp collisions than in heavy ion collisions, e.g. 
the freeze-out chemical potential follows a different pattern, due to the lower stopping power in pp 
collisions. 

The relation between the p/p ratio and Hb can be shown easily within the statistical concept 



using the Boltzmann statistics Ref. []33[]. In the model calculation, the appropriate statistics and 
also feed down from strong decays are taken into account. The density of particle i is then given by 



di 



2n 2 V T J 



with Nb and Ns being the baryon and strangeness quantum numbers of particle i. 
This leads to a p/p ratio of (excluding feed-down from heavier resonances): 



(4.3) 



(4.4) 



The ratio of strange antibaryons/ baryons is then given by 



% =e -(2ii B -N sl i s )/T 

n B 



(4.5) 
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Figure 5: Variation of the baryon chemical potential fig as a function of \J~s. The dashed line describes 
heavy ion collisions as in Ref. [ p2| | while the solid line is the new parametrization for pp collisions. 

As Us is always smaller than /!#, the ratios appear ordered with the strangeness quantum 
number, i.e. the higher N$, the smaller the difference between antibaryon and baryon. This trend 
is shown in Figs. ^ and ^ comparing the results from the model with experimental data. The 
agreement between the model results and the data is very good. 

5. Production of nuclei, antinuclei, hypernuclei and antihypernuclei 



5.1 Comparison to data from RHIC 

The production of light nuclei including hypertritons (^H) and antihypertritons (^H) was re- 
cently observed by the STAR collaboration JTr^]. The abundances of such light nuclei and antinuclei 
follows a consistent pattern in the thermal model. The temperature remains the same as before but 
an extra factor of /ig is picked up each time the baryon number is increased. Each proton or neutron 
thus simply adds a factor of /!# to the Boltzmann factor. The production of nuclear fragments is 
therefore very sensitive to the precise value of the baryon chemical potential and could thus lead to 
a precise determination of /!#. 

The ratios within the statistical approach using the grand-canonical formalism can be easily 



written, based on Eq. ( |4.3| ). Deuterium has an additional neutron and the antideuterium to deu- 
terium ratio is given by the square of the antiproton to prton ratio: 



-(*««)/r 



"d 



(5.1) 
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Figure 6: Antibaryon to baryon ratios at the SPS according to strangeness content. Circles refer to p- 
p collisions, squares to heavy ion collisions. 




Figure 7: Antibaryon to baryon ratios at STAR according to strangeness content. Circles refer to p-p colli- 
sions, squares to heavy ion collisions. 
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Helium 3 has 3 nucleons and the corresponding anti-Helium 3 to helium 3 ratio is given by: 

(5.2) 



^ = e -(6na)/T 



If the nucleus carries strangeness this leads to an extra factor of jj, s 



m-. 



JlI =e -(^B-2ns)/T 



«3, 



(5.3) 



In mixed ratios the different degeneracy factors are also taken into account, e.g. 6 for 3 A H and 2 for 



I*. 



m i 



3 e -(6/te-Ms)/r 



m 



(5.4) 



He 



In the model like in the data the He 3 and He 3 yields have been corrected for the part coming from 
hypertriton and antihypertriton decays assuming a decay branch ratio for the decay of 25 %. 

5.2 Predictions for RHIC and LHC 

In Fig. H we compare p-p and heavy ion collisions at ^/s = 200 GeV. The difference between 
the two colliding systems and the effect of canonical suppression is seen in p-p collisions. 

In Fig. 7 a comparison is shown of the various antiparticle/partcle ratios for two different beam 
energies. 



The expectations for the LHC are shown in Fig. 10 
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Finally the predictions of the thermal model for ratios of anti-nuclear to nuclear fragments are 
shown in Fig. [0] This figure includes comparisons for strange nuclear fragments where a clear 
picture emerges (again) between strange and non-strange fragments. 
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Figure 11: The ratio of the yield for examples of different masses. 

6. Discussion and Summary 

In the present paper we have made a general comparison of thermal parameters in p-p and 
heavy ion collisions. We have determined the energy dependence of the baryon chemical potential 
/ifiin p-p collisions. This was used to establish a hierarchy of antibaryon to baryon ratios includ- 
ing strange and multi-strange baryons. This was then used to compare nuclear and anti-nuclear 
fragments in p-p and heavy ion collisions. Predictions have been presented for these ratios at LHC 
energies. 
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